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ABSTRACT 

Aims. Our aims are twofold. First we aim to evaluate the robustness and accuracy of stellar parameters and detailed elemental 
abundances that can be derived from high-resolution spectroscopic observations of microlensed dwarf and subgiant stars. We then 
aim to use microlensed dwarf and subgiant stars to investigate the abundance structure and chemical evolution of the Milky Way 
Bulge. Contrary to the cool giant stars, with their extremely crowded spectra, the dwarf stars are hotter, their spectra are cleaner, and 
the elemental abundances of the atmospheres of dwarf and subgiant stars are largely untouched by the internal nuclear processes of 
the star. 

Methods. We present a detailed elemental abundance analysis of OGLE-2008-BLG-209S, the source star of a new microlensing event 
towards the Bulge, for which we obtained a high-resolution spectrum with the MIKE spectrograph on the Magellan Clay telescope. 
We have performed four different analyses of OGLE-2008-BLG-209S. One method is identical to the one used for a large comparison 
sample of F and G dwarf stars, mainly thin and thick disc stars, in the Solar neighbourhood. We have also re-analysed three previous 
microlensed dwarf stars OGLE-2006-BLG-265S, MOA-2006-BLG-099S, and OGLE-2007-BLG-349S with the same method. This 
homogeneous data set, although small, enables a direct comparison between the different stellar populations. 

Results. We find that OGLE-2008-BLG-209S is a subgiant star that has a metallicity of [Fe/H] « -0.33. It possesses [or/Fe] 
enhancements similar to what is found for Bulge giant stars at the same metallicity, and what also is found for nearby thick disc stars 
at the same metallicity. In contrast, the previous three microlensing dwarf stars have very high metallicities, [Fe/H] > +0.4, and more 
solar-like abundance ratios, i.e. [ff/Fe] « 0. The decrease in the [ff/Fe] ratio with [Fe/H] is the typical signature of enrichment from 
low and intermediate mass stars. We furthermore find that the results for the four Bulge stars, in combination with results from studies 
of giant stars in the Bulge, seem to favour a secular formation scenario for the Bulge. 

Key words, gravitational lensing - Galaxy: bulge Galaxy: formation - Galaxy: evolution - Stars: abundances - Stars: fundamental 
parameters 



1. Introduction 

How spiral galaxies form and acquire their different stellar pop- 
ulations is largely an unsolved problem. However, since the 
mid-1980s the theory of hierarchical structu re formation in a 
Lamb da Cold Dark Matter Universe (ACDM, Blumentha l et all 
1984) has emerged as a foundation to understand the proper- 
ties and evolution of galaxies. In ACDM cosmological simu- 



* Based on data collected with the 6.5 m Magellan Clay telescope at 
the Las Campanas Observatory, Chile. 



lations, galaxies such as the Milky Way are built from inde- 
pendent smaller systems and fragments/debris from other stel- 
lar systems over a time period spanning a few billion years (e.g. 
iGovernato et al .1120071 ; iRead et al.ll2008l) . A fundamental predic- 
tion from such models would be that bulges formed in merg- 
ers. However, at the same time there is growing observational 
evidence that the bulges of man y distant galaxies fo rmed by in- 
ternal dynamical processes (e.g jGenzel et a l. 2008). During an 
early turbulent phase of a galaxy, disc material (gas and stars) 
is gravitationally driven into the central regions of the galaxy, 
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building up an exponential c omponent. This scenario is re ferred 
to as secular evolution (e.g., Kor mendv & Kennic utt 2004). 

The Milky Way bulge (hereafter referred to as the Bulge, 
with capital B) has a rather broad metallicity distribution (e.g., 
IZoccali et al.l 12008) and its metal -rich stars a nd globular clus- 
ters are as old as the G alactic halo (e.g., iRosenberg et alj 
1 19991; iMarin-Franch et al.|[20Q8b . This points to a very intense 
star formation rate in the early history of the Galaxy. The 
Bulge thus represents an important link for our understand- 
ing of galactic bulges in general, and is int egral to the ques- 
tion of galaxy formation and evolution (e.g.. IWvse et al.| [l997: 
iKormendv & Kennicutl2004l) . 

Despite the fact that detailed abundance studies can pro- 
vide crucial information on the formatio n and chemical en rich- 
ment history for a stellar population (e.g. lMcWilliam|[l997l) . the 
Bulge has long been the least studied stellar component of our 
Galaxy. This is due to the inherent difficulties in studying stars 
in the Bulge as they are distant and suffer from a high degree 
of interstellar extinction in th e Galactic plane. However , fol- 
lowing the pioneering study by iMc William & Rich! {l994), and 
with the advent of 8-10 meter class telescopes, substantial in- 
sight into the stellar populations and the chemical history of the 
Bulge has been gained from high-r esolution spectroscopic stud- 
ies of bright K and M giant stars dFulbright et alj|2006l 120071: 



Cunha & Smith 2006; Cunha et al. 2007, 2008; Rich & Origlia 
20051: i Rich et"al.ll2007j;lLecureur et al.ll2007HZoccali et al.l l2003. 
2008; lMelendez eta l. 2008; Rvde et al.ll2009l) . 



When studying the chemical history of a stellar population 
through detailed elemental abundances, one relies on the as- 
sumption that the chemical composition of the stars is a true 
measure of the elemental abundances present in the gas from 
which they formed. The expected lifetimes for F and G dwarf 
stars on the main sequence, burning hydrogen to helium in their 
centres, are similar to, or possibly even longer than the current 
age of the Galaxy. For instance, a sol ar type star will spend 
~ 11 Gyr on the main sequence (e.g. ISackmann et al.l 1 1993b 
during which their atmospher es are un touched by internal nu- 
clear processes of the star (e.g. Iben 1991). At later evolutionary 
stages, when the stars reach the red giant branch, various inter- 
nal physical processes can erase the abundance signatures from 
the stellar atmosphere. This is the case for C, N, and Li in all 
giant stars , and for O, Na, Mg , Al in giant stars in globular clus- 
ters (e.g. [Gratton et al. 2004). For dwarf stars it is essentially 
only the rare lig ht elements Li, Be and B that are depleted (e.g., 
Boesgaard 2005), and as studies of chemical evolution mostly 
focus on heavier elements, such as the a-elements, this is one of 
the main reasons why observational studies of the chemical evo- 
lution of the Milky Wa y are largely based on dwarf stars (e.g. 
Edvardss on et aljfl993l) . Also, because the spectra of the intrin- 
sically luminous metal-rich giant stars have a very rich and often 
hard to identify fauna of molecular bands, great effort has to be 
made to identify un -blended and weak spectral lines (see, e.g., 
iFulbright et al.ll2007l) . Furthermore, for species such as, e.g., Mg 
and Na, that only have a few usable lines that usually also are 
strong (> 1000 mA) in metal-rich dwarf stars, the lines become 
very strong in giant stars. Elemental abundance uncertainties due 
to NLTE effects are also likely to be larger in giant stars than in 
dwarf stars (e.g., lAspiundll2005l) . Hence, there is reason to be- 
lieve that the underlying assumption that the giant stars accu- 
rately trace the chemical evolution of a stellar population could 
be erroneous and should therefore be rigourously tested. 

The concept of using microlensing to obta in spectra of 
dwarf stars in the Bulge was first demonstrated by M inniti et alj 
(1998) who used Keck I as a "15 meter telescope" to obtain 
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Fig. 1. Dwarf stars in the Bulge for which high-resolution spectra 
have been obtained. Filled circles mark OGLE-20 06-BLG-265 S 
and M OA- 2006-BLG-099S alre ady published bv lJohnson etail 
d2007l) and lJohnson et al.l d2008l), respectively , and OGLE-2007- 
BLG-349S published by ICohen etal] (12008b . The "star" marks 
OGLE-2008-BLG-209S being analysed here for the first time, 
and the open circles mark OGLE-2007-BLG-514S, MOA-2008- 
BLG-310S, and MOA-2008-BLG-31 IS that currently are being 
analysed elsewhere. The curved lines show the outline of the 
Bulge based on observations with the COBE satellite. The open 
square marks the position of Baade's window (BW), and the 
Galactic centre (GC) is also marked. 



a s pectrum of th e mod erately magnified (~ 1 mag) 97-BLG- 
45. ICavallo et al.l (12003) then presented an analy sis of six mi- 
crolen sed stars (including the star observed by iMinniti et al.l 
119981) : two cool giant stars, one subgiant star, t wo solar ana- 
logue s, and one of uncertain type (the one from IMinniti et al] 
179981) . 

After the first studies by IMinniti et aT] (1 19981) and 
ICavallo et al.1 (120031). three microlensed Bulge dwarf stars hav e 
been observed dJohnson etani2007l l2008t ICohen et al.l 12008). 
The results from these studies provide some surprising results, 
contradicting the results based on giant stars. First, it is clear 
that all three are un usually metal-rich. The star analysed by 
iJohnson et all (120071) . at a metallicity of [Fe/H] = +0.56, is the 
most metal-rich star known to us. As microlensing events have 
no bias with regard to the metallicity of the source star, analysing 
enough events will provide an unbiased metallicity distribution 
(MDF) of the Bulge. The three dwarf events studied in detail 
so far point to a much more metal-rich MDF compared to the 
one derived from giant stars. Secondly, the stars have elemen- 
tal abundance ratios that in general are very similar to what is 
found in metal-rich thin disc stars in the Solar neighbourhood, 
i.e. they do not show the high [a/Fe] ratios that is found in the 
Bulge giant stars. Although the dwarf sample is still very small, 
these findings hint that the assumption that giant stars give us a 
complete picture of the Bulge may not be correct. 

In this paper we present a detailed abundance analysis of 
OGLE-2008-BLG-209S, the source star of a microlensing event 
towards the Bulge. We have performed four different analyses of 
the star, enabling a comparison of the robustness of the derived 
stellar parameters and elemental abundances. The methods we 
use have all been applied either in previous studies of of lensed 
Bulge dwarf stars or in large studies of nearby dwarf stars. One 
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method utilises the spectral linelist, atomic data, model stellar 
atmospheres, and method to find the stellar parameters that are 
all very sim ilar to what is used in the studies of nearby F and G 
dwarf stars dBensbv et al.ll2003ll2005L and Bensby et al. 2009, in 
prep.). This enables direct comparisons between the Bulge stars 
and the local disc stars. The other two methods are similar to the 
ones used for the analyses of OGLE-2006-BLG -265S, MOA- 
2006-BLG-099S, and OGLE-2 007-BLG-349S in dJohnson et alj 
120071, l2008HCohen et alj2008l) . Spectra for these three other mi- 
crolensed Bulge dwarf stars were made available for this study 
and we have re-analysed them wit h the method akin to the one 
used in iBensbv et ail J2003L [2005b and Bensby et al. (2009, in 
prep.). 

2. Observations and data reduction 

In Fig. [1] we show the positions of the microlensed dwarf 
stars in the Bulge that have been observed with high-resolution 
spectrographs so far. In this figure, three additional events are 
marked. These three are OGLE-2007-BLG-514S, MOA-2008- 
BLG-310S, and MOA-2008-BLG-311S. They are currently be- 
ing analysed and will be published elsewhere (Epstein et al., in 
prep.; Cohen et al., in prep.). 

2.1. Nomenclature 

First we want to clarify the nomenclature used to describe mi- 
crolensing events. The actual microlensing event is given a name 
depending on the project which discovered it (e.g., MOA or 
OGLE), year it took place, the direction in which the event oc- 
curred (in our case BLG for the Bulge), and a running number, 
e.g., OGLE-2008-BLG-209. When referring to the source star of 
the microlensing event an "S" is added to the event name, i.e., 
OGLE-2008-BLG-209S. Other extensions are given for the lens 
and for any planets orbiting the lens. 

2.2. OGLE-2008-BLG-209S 

In May 2008 the OGLE early warning systerrfl dUdalskill200"3h 
put out an alert that a new microlensing event, OGLE-2008- 
BLG-209, was developing and would reach maximum bright- 
ness on HJD2454606.097. The faintness of the source star be- 
fore the microlensing event made it likely for the source star to 
be a dwarf star in the Bulge, located close to Baade's window 
at (/, b) = (1.4, -4.0) deg (see Fig.Q]). Fortunately we had ob- 
serving time on Magellan II telescope at this time allowing us to 
obtain a high-r esolution spectrum of the object with the MIKE 
spectrograph dBernstein et al.ll2003l) . At the time of observation 
on May 20, 2008 (HJD2454606.838) the event had just passed 
maximum brightness and had a magnitude of / ~ 14.5, about 
20 times brighter (approximately 3 magnitudes) than before be- 
ing microlensed (see Fig. |2j. Three 1800 s exposures were ob- 
tained using a slit width of 0.7 arcsec, resulting in a spectrum, 
with continuous wavelength coverage from 3200 A to 10000 A, 
recorded on two CCDs. With a 0.7 arcsec slit the resolving power 
is R ~ 55 000 on the blue CCD and R ~ 47 000 on the red CCD. 
The data were reduced with the Carnegie Observatories MIKE 
Python pipeline (D. Kelson, private communication) and the fi- 
nal spectrum used in the analysis consists of the three individual 
spectra co-added and then divided by the blaze function. Typical 
signal-to-noise ratios are S /N ~ 20 per pixel in the spectrum 



1 http : //ogle . astrouw . edu . pl/ogle3/ews/ ews . html 
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Fig. 2. OGLE photometry of the microlensing event. Maximum 
brightness was estimated to occur on HJD2454606.097 (+0.004) 
i.e., UT2008-05-20.06. The solid fine is a theoretical fit to the 
microlensing event. The vertical dotted line indicates when ob- 
servations (start of first exposure) were carried out with MIKE 
(HJD2454606.838). 



from the blue CCD and S /N ~ 30 per pixel in the spectrum 
from the red CCD (as measured in the continuum of regions free 
of spectral lines). 

2.3. OGLE-2006-BLG-265S, MOA-2006-BLG-099S, and 
OGLE-2007-BLG-349S 

iJohnson et alj d2007l) presented a detailed abundance analysis of 
OGLE-2006-BLG-265S. At the time of observation it was mag- 
nified by a factor of ~ 145 and, by using the HIRES spectro- 
graph on the Keck telescope, they obtained a spectrum of high- 
resolution (R ~ 45 000) and relatively high signal-to-noise ratio 
(S IN ~ 60 per resolution element). 

A year later. I Johnson et all (1200 8) presented a detailed abun- 
dance analysis of MOA-2006-BLG-099S, another microlensed 
Bulge dwarf star. This object was at the time of observation mag- 
nified by a factor of ~ 110, and a spectrum was obtained with 
the MIKE spectrograph on the Magellan II telescope. Compared 
to OGLE-2006-BLG-265S, the spectrum of MOA-2006-BLG- 
099S is of lower quality with a resolution of only R ~ 19 000 
red-ward of 4800 A and S /N ~ 30 per pixel. 

OGLE-2007-BLG-349S, the third dwarf reanalysed here, 
was observed on September 5, 2007 at the Keck Observatory. 
Three consecutive spectra, each 1350 s in length, were obtained 
with HIRES-R at the Keck I telescope in a configuration with 
coverage from 3900 to 8350 A, with small gaps between the 
orders beyond 6650 A. The slit was 0.86 arcsec wide, giving a 
spectral resolution of R ~ 48,000. The magnifications at the 
time of the recording of the three spectra were 350, 390, and 
450. The signal-to-noise ratio of the resulting summed spectrum 
is S IN ~ 30 per spectral resolution element for wavelengths be- 
low 5000 A and S/N -90 for wavelengths above 5500 A. 



4 



T. Bensby et al.: Chemical evolution of the Galactic bulge as traced by dwarf and subgiant stars 



— , 150 



100 



50 



1 1 1 1 1 1 1 1 1 1 


1 1 1 | 1 £ 


- OGLE-2006- 




- BLG-265S 












/ 


-J=-2% " 


r B 

•i i i i i i i i i i 


(9%) - 
i i i 1 i i 



50 100 150 
^Bensby [ mA ] 



,— , 150 



100 



50 



1 1 1 1 1 1 1 1 1 1 
- OGLE-2008- 




- BLG-209S 




\/ 


B-J=+2% I 


1 1 1 1 1 1 1 1 1 


(9%) - 
i i i i 1 i i 



50 100 



150 



Bensby 



[mA] 



^ 150 



100 



e 



c 

o 



50 



i i i i I i i i i I i i i i I i i 

- MOA-2006- • • .•' 

- BLG-099S , • ••• 

/ 

B-J=-4% 
(14%) 

•i i i i I i i i i I i i i i I i i 



50 100 150 
^Bensby [ mA ] 



150 



100 



50 



i i i i I i i i i I i i i i I i 
OGLE-2008- 
BLG-209S 



J-C=+3% 
* (8%) 
i i i i I i i i i I i i i i I i i 



50 100 150 
^Johnson [ mA l 



J3 

O 



c 

o 



150 



100 



50 



1 1 1 1 1 1 1 1 1 1 1 


, , , | Mj 


- OCLE-2007- 




- BLG-349S 


m 


/ 


-C=-l% " 


'i i i i 1 i i i i 1 i 


(5%) - 
i i i 1 i i 



50 100 150 
^Bensby t mA ] 



150 



100 



50 



i i i i I i i i i I i i i i I i i 
OGLE-2008- 
BLG-209S 



B-C=+5% 

*r (n%) 
i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



50 100 150 
^Bensby t mA ] 



Fig. 3. Comparison of equivalent widths measured independently by Bensby (B), Cohen (C) and Johnson (J) for the four stars that 
have been analysed. Differences and standard deviations (in parentheses) are indicated in the figure. The dotted line in each plot 
indicates the one-to-one relation. 



Further details about obser vations and data reduction s for 
these three stars can be found in iJohnson et ail d2007ll2008l) and 
ICohen et al.1 d2008l) . 



2.4. Radial velocities 

For OGLE-2008-BLG-209S we measure a velocity shift of 
-190.1 ± 0.4 kms -1 in the spectrum from the blue CCD and 
-189.0 ±0.1 kms" 1 in the spectrum from the red CCD. As 
the red spectrum is of higher quality we use this one only for 
the radial velocity. Adding the heliocentric correction, which 
at the time of observation was + 15.4 km s~', the heliocen- 
tric radial velocity of OGLE-2008-BLG-209S is -173.6 kms 1 . 
The three stars, OGLE-2006-BLG-265S, MOA-2006-BLG- 
099S, and OGLE-2007-BLG-349S ha ve radial velocities of +99 
-154, and +113 k m s _1 , respectively ( I Johnson et al.l2007l l2008t 
ICohen et al.ll2008l) . 



3. Abundance analysis 

Stellar parameters and elemental abundances will be determined 
independently by three of us (T.B., J. J., and J.C.) using different 
approaches. The different methods will use their own linelists, 
atomic data, and model stellar atmospheres as used in the various 
previous studies of Bulge as well as local dwarf stars. As equiv- 
alent widths are the fundamental ingredient for all four methods, 
we will start by checking the agreement between the different 
measurements for lines that are in common. 



3.1. Linelists and equivalent width comparisons 

The analysis by Bensby (methods 1 a nd 2 below) uses the sam e 
spectral line list and atomic data as in ((Bensbv et al.l2003ll2005l) . 
However, it has now been expanded by another 50 Fe i lines (see 
Bensby et al. 2009, in prep.; and Table J2|. Equivalent widths 
were measured by hand using the IRAHj task splot. Gaussian 
line profiles were fitted to the observed lines and in special cases 
of strong Mg, Ca, Si, and Ba lines, we used Voigt profiles instead 
to better account for the extended wing profiles of these lines. 

The analysis by Johnson (metho d 3 below) use s a linelist tha t 
for most elements was taken from Bensb vet alj (120031 [2004). 
The equivalent widths were measured using spectrjUIc. 
Sneden, 2007, private communication). For Ba, synthetic spectra 
were compared to the observed spectra to determine the abun- 
dance. The effect of hyperfine splitting (HFS) in the Ba lines 
was included and the Ba HFS constants and log g f- v alues were 
taken from the sources listed in I Johnson et alj d2006l) . 

The analysis by Cohen (method 4 below) uses the linelist 
given in Cohe n etaT] d2008l) . Equivalent widths were measured 
using an automatic Gaussian fitting routine, after which stronger 
lines were checked by hand in order to make sure that damping 
wings were picked up when appropriate. Elements with only a 
few detected lines were also checked by hand. 

Figure [3] shows comparisons of the equivalent widths as 
measured by Bensby, Johnson, and Cohen. There is generally 
a good agreement. For OGLE-2008-BLG-209S Bensby's mea- 

2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under co-operative agreement with 
the National Science Foundation. 

3 http : //verdi . as . utexas . edu/ spectre . html 
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Table 1. Comparison of equivalent widths measured by Bensby, Johnson and Cohen. For each element we give the total number of 
lines that each of the different linelists contain, the average differences (in percent), the 1-sigma spreads, and the number of lines in 
common (in parentheses). 





OGLE-2008-BLG-209S 
Ben John Coh Ben-John Ben-Coh John-Coh 

N N N %±cr(N) % ±cr (N) % ±cr (N) 


MOA-2006-BLG-099S 
Ben John Ben-John 

N N % ±cr (AO 


OGLE-2006-BLG-265S 
Ben John Ben-John 

N N % ±cr (A?) 


OGLE-2007-BLG-349S 
Ben Coh Ben-Coh 

N N % ±cr (N) 


Total 
Fei 
Fe ii 


284 146 109 +2 ±9 (135) +5± 11 (55) +3 ± 8 (32) 
145 72 101 -1 ± 10 (61) +4±11(49) +3 ±8 (31) 
19 2 8 +4 ±8 (2) +12 ±12 (6) +9 ±0(1) 


251 105 -4 ±14(76) 
122 35 -5 ±11 (19) 
14 4 -3 ± 9 (3) 


188 219 -2 ±8 (172) 
91 98 -2 ±8 (79) 
21 9 -3 ±14 (8) 


181 249 -1 ± 5 (1 17) 
91 135 -0 ± 3 (59) 
14 11 -2 ± 8 (11) 



surements are on average 5 % larger than Cohen's (55 lines in 
common), and 2% larger than Johnson's (135 lines in com- 
mon), which in turn are 3 % larger than Cohen's (34 lines in 
common). For OGLE-2006-BLG-265S the Bensby measure- 
ments are just 2 % smaller (172 lin es in common) than the mea- 
surements by iJohnson et at] d2007l) . and for MOA-2006-BLG- 
099S the Bensby measurem ents are 4% s malle r (76 lines in 
common) than the ones by IJohnson et al.l d2008l) . For OGLE- 
2007-BLG-349S, which has the highest SNR spectrum of the 
four microlensed stars, the agreement is even better, the Bensby 
measuremen t s are only 1 % smaller than the measurements by 
ICohen et al.l {2008) (117 lines in common). These differences 
relate to all equivalent widths from all elements and ions mea- 
sured. By comparing the equivalent widths of the Fei and 
Fen lines only, we see that there are significant larger aver- 
age differences for the Fe n lines than for the bulk of lines (see 
Table [Q. For instance, for OGLE-2008-BLG-209S the differ- 
ences in equivalent widths between Bensby's and Cohen's mea- 
surements of Fe ii lines reach 12%. How these differences affect 
the stellar parameters is further investigated in Sect. 13.61 

Figure [4] shows a comparison of oscillator strengths (log gf- 
values) for Fei and Fen lines. Since the Johnson studies 
largely have a dopte d the linelists and atomic data as given in 
Bens by et al.l d2003l) we only make comparisons between the 
Bensby and Cohen linelists. The figure shows the differences 
in the log gr/-values as a function of wavelength, as a function 
of the lower excitation potential of the line, and as a function 
of the equivalent width. From these plots we see that there are 
differences present, however with no obvious trends with either 
excitation potential, line strength, or wavelength. Average differ- 
ences are -0.03 ± 0.07 dex for Fe i lines and -0.05 ± 0.09 dex for 
Fe n lines, with the values from Cohen linelist being the larger 
ones. This means that for a given combinati on of stellar parame- 
ters, the average absolute iron abundances bv lCohen et all (120081) 
will come out 0.03 dex lower if based on Fe i lines and 0.05 dex 
lower if based on Fe ii lines, if the same equivalent widths and 
model atmospheres are used. 

As the tuning of the stellar parameters are based on absolute 
abundances from Fe i and Fe n lines (i.e. not normalised to the 
solar abundances) these differences may have an effect on how 
the stellar parameters come out from the different methods. This 
will be investigated further in Sect. 13.61 Note that it is essentially 
only for the Fe lines that good accuracy of the log ^/-values are 
important. Abundances based on lines from other elements will 
be normalised to those of the Sun, cancelling, to first order, out 
any uncertainties in the log ^/-values. 

Also, the two Fei lines at 6726 A and 7941 A clearly show 
larger deviations than the other lines (see Fig. @}. In Table [2] 
where we give abundances for all lines we see that the solar 
abundance from the /16726 line is too low, and from the A1941 
line it is too high. The normalised abundances from these two 
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Fig. 4. Comparison of Fe i and Fe n oscillat or strengths (log q f) 
used in methods 1 and 2 (values taken from lBensbv et al.l f2003 



and Bens by et al. 2009 , in prep.) and method 4 (values taken 
from ICohen et al.l 12008). The plots show Bensby's values mi- 
nus Cohen's values as a function of line strength (top panel), 
in this case Cohen's measurements for OGLE-2007-BLG-349S; 
lower excitation potential (middle panel); and wavelength (bot- 
tom panel). Open and filled circles mark Fei and Fen lines, 
respectively. On average the Cohen values are 0.03 + 0.07 dex 
larger for the Fei lines and and 0.05 + 0.09 dex larger for the 
Fe ii lines (i.e. giving lower abundances). 



lines for the four stars agree however well with the normalised 
abundances from the other Fe i lines. So, it seems that the log gf 
values for these two lines could be significantly off. However, as 
disregarding the two lines have no effects on the stellar param- 
eters nor the final abundances, we will, for the moment, keep 
them as they are. 

The Bensby linelist, with measured equivalent widths and 
calculated elemental abundances for OGLE-2008-BLG-209S, 
OGLE-2006-BLG-265S, MOA-2006-BLG-099S, and OGLE- 
2007-BLG-349S is given in Tabled 

3.2. Method 1: Bensby analysis with gravity from ionisation 
balance 

This is a standard method based one-dimensional plane- 
parallel LTE model stellar atmospheres that were calcu- 
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Fig. 5. Absolute Fe abundances versus reduced line strength and versus lower excitation potential of the line. Empty circles are 
abundances from Fe i lines and filled circles from Fe n lines. The dashed lines are the linear regression to the Fe i abundances. 



lated with the Uppsala MARCS code dGustafsson et alj|1975t 
lEdvardsson et alj|1993t lAsplund et a"T]|1997l) . Elemental abun- 
dances were calculated with the associated suite of programs 
for abundance analysis and line synthesis (EQWIDTH and 
SPECTRUM, Edvardsson, Eriksson, & Gustafsson, 2000, pri- 
vate comm.). The resulting stellar parameters and elemental 
abundances are all based on the equivalent widths measured by 
T. Bensby. The choice of these MARCS model atmospheres en- 
sures that the results based on this method are fully compatible 
with our studies of ne arby dwarf stars bel onging to the Galactic 
thin and thick discs (Bensb v~et al.ll2003l 12005, and Bensby et 
al. 2009, in prep.). To find the stellar parameters we use a grid of 
approximately 12 000 MARCS model stellar atmospheres span- 
ning metallicities between -2.2 < [Fe/H] < +0.7 in steps of 
0.1 dex; surface gravities between 3.2 < \ogg < 4.8 in steps 
of 0. 1 dex, and temperatures between 4500 < T e g < 7000 K 
in steps of 100 K (see Bensby et al. 2009, in prep.). The mod- 
els have chemical com positions sca l ed rela tive to the standard 
solar abundances of Asplu ndet alj d2005l) . However, to bet- 
ter reflect the actual composition of the stars, starting at solar 
metallicity, the models have a-enhancements that increase with 
decreasing metallicity, reaching [a/Fe] = 0.40 at [Fe/H] = 
-1.0. Below [Fe/H] = -1.0 the (^-enhancement is constant at 
[a/Fe] = 0.40, and above the solar metallicity, [Fe/H] = 0, there 
is no a-enhancement. When calculating ab undances the broad - 
ening by collisions were taken fr o mlAnstee & O'Maral (Tl995); 
IBarklem & O'Maral d 19971 Il998h : iBarklem et alJ (119981 l200Q>. 
For lines not includ ed in these works the classical Unsold broad- 
ening was used (see Bensby et al. 2003b . 

The stellar parameters are determined using the abundances 
from Fei and Fen lines. The basic concepts of the method are 



as follows: (1) The effective temperature (T s ff) is determined by 
requiring zero slope in the diagram where abundances from the 
Fei lines are plotted versus the lower excitation energy of the 
line (x\), i.e. excitation equilibrium; (2) The microturbulence 
parameter (£ t ) is determined by requiring zero slope in the di- 
agram where abundances from Fe i lines are plotted versus the 
measured line strength (log(W//l)); (3) The surface gravity is de- 
termined by requiring that the average abundance from the Fe i 
lines equals the average abundance from Fe n lines, i.e. ionisa- 
tion equilibrium. This method was used because we do not know 
the apparent magnitude nor the distance of the star. All of these 
balances are coupled and it is not possible to first tune one pa- 
rameter and then another. All parameters have to be found simul- 
taneously, and after extensive testing we find that an acceptable 
solution can only be found for one set of parameters for a given 
star. Final diagnostic plots for and ^ t are shown in Fig. [5] The 
ionisation balance was deemed good when the average output 
abundances from the Fe i and Fe n lines agreed within 0.01 dex. 
Also, to avoid uncertainties in the stellar parameters arising due 
to saturation effects in strong lines, only Fe i and Fe n lines with 
measured equivalent widths smaller than 90 mA were used in the 
determination of the stellar parameters. 

There are indications that abundances from some Fe i lines 
may be sensitive to departures from LTE which could invalidate 
our assumption of ionisation balance when determining the sur- 
face gravity. However, the predicted correcti on for a dwarf star 
at fFe /H] = -0.5 is s mall, around . 05 dex ( Theveni n & Idiart 
1999). Furthermore, in iBensbv etaTJ d2003l 120051) and Bensby 
et al. (2009, in prep.) we have analysed a total of ~ 700 nearby 
F and G dwarf stars in the Solar neighbou rhood, all hav i ng ac- 
curate distances from the new reduction by Ivan Leeuwenl (2007) 
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of the parallaxes from the Hipparcos satellite. The fact that for 
these ~ 700 F and G dwarf stars we do not require ionisation 
equilibrium, but anyway find a very good agreement between 
abundances from the two ionisation stages of iron ([Fe i/Fe n] = 
0.01 + 0.07), means that ionisation balance should be a valid 
option to use for F and G dwarf stars when the distance and/or 
the apparent magnitude is unknown. Hence, Met hod 1 and the 
method for nearby dwarf stars as presented in Bensb v~et al.l 
d2003l l2005h and Bensby et al. (2009, in prep.), where the dis- 
tances and magnitudes of the stars are known, should produce 
results that are fully compatible. 



3.3. Method 2: Bensby analysis with gravity based on 
microlensing assumptions 

This method is similar to method 1 and only differs in the 
way the surface gravity is estimated. Instead of using ionisa- 
tion balance the surface gravity will be determined using an ab- 
solute magnitude estimated from microlensing techniques. As 
the linelists, equivalent width measurements, model stellar at- 
mospheres, are the same as in method 1, we only describe how 
the absolute magnitude is estimated. 

De-reddened colours and magnitudes of the source can be es- 
timate d using standard microlensing techniques (e.g. lYoo et all 
|2004|) . The method for determining the colour does not make any 
assumption about the absolute reddening, nor about the ratio of 
selective to total extinction. It only assumes that the reddening 
to the microlensed source is the same as the reddening to the red 
clump, and that the red clump has (V - 7)o = 1 .00, the same as 
the local Hipparcos red clump. However, in principle, the Bulge 
red clump could be different due to different age and composi- 
tion of the Bulge as compared to the Solar neighbourhood stars. 
From the sp ectroscopic temperatures of the three previous B ulge 
dwarf stars dJohnson et al.ll2007l 120081: ICohen et al]|2008l) it is 
found that the mean Bulge clump has (V-I)o = 1.05. This value 
is in good agreement with a completely independent estimate 
based on taking into account the age and differences in composi- 
tion between the Bulge and the Hipparcos stars (David Bennett, 
private communication). (V - I)o = 1.05 for the red clump will 
be the value used in this study. The absolute de-reddened colour 
is then derived from the colour offset between the microlensing 
source and the clump in the instrumental colour-magnitude dia- 
gram (CMD). V and I measurements give a colour estimate of 
(V - 7)o = 0.73 for OGLE-2008=BLG-209S. 

The de-reddened apparent magnitude of the source is de- 
termined in a similar way, using an assumed de-reddened ap- 
parent magnitude of the red clump (7o = 14.32) and the offset 
between the apparent magnitudes of the source and the clump 
on the instrumental CMD. The instrumental magnitude of the 
source is determined from the microlensing model. This yields 
7 = 17.05 for OGLE-2008-BLG-209S. The absolute magni- 
tude of the source is then determined from assuming that the 
red clump stars are at the same distance as the source. For 
OGLE-2008-BLG-209S the estimate for the apparent magnitude 
is Vrj = 17.68, and the absolute magnitude is then My = 3.26 
(calculated assuming a distance of 8 kpc). 

Using this absolute magnitude the stellar parameters are then 
found in the same way as for method 1, except that ionisation 
equilibrium is not required. 



3.4. Method 3: Analysis by Johnson 

TurboSpectrum ( I Alvarez &Plezl [[998), a 1 -dimensional LTE 
code, was used to derive elemental abundances. The input 
model stellar atmospheres were interp olatecfl in a grid of th e 
new MARCS 2008 model atmospheres (Gustafs son et al.ll2 008). 
These models have chemical compositions scaled r elative to the 
standard solar abundances in lGrevesse et al.l d2007t) . but with a- 
enhancements for sub-solar metallicities that are the same as 
used for the MARCS 1997 models in methods 1 and 2 above (see 
Sect. |3.2| >. When calculating elemental abundances the treatment 
of damping from Barkler rTet al.l d2000l) was used. 

Effective temperatures and the microturbulence velocities 
were determined using standard techniques in a similar way to 
Method 1 , including log g that was determined by demanding 
ionisation equilibrium. The final parameters adopted for OGLE- 
2008-BLG-209S were r eff = 5250 K, \ogg = 3.5, metallicity 
[Fe/H] = -0.40, and£ = 1.4 kms" 1 . 

Further details of the method can be found in [Johnson et al.l 
( 20071 [2008h . It should be noted that an attempt to measure the 
effective temperature fro m the Balmer Ha line, as was done for 
MOA-2006-BLG-099S in lJohnson etail (120081) . was made, but 
the cool temperature of OGLE-2008-BLG-209S meant that no 
accurate T e $ could be derived this way. 

3.5. Method 4: Analysis by Cohen 

This method is also based on standard LTE assumptions, using 
a grid of solar sca l ed AT LAS9 model stellar atmospheres from 
ICastelli & Kuruczld2003l) . In the analysis for OGLE-2008-BLG- 
209S, only lines with wavelengths longer than 5400 A were se- 
lected due to the low S /N in the crowded blue parts of the spec- 
trum. Also, the stronger lines in these high gravity dwarf stars 
with their extended damping wings, which are difficult to mea- 
sure in crowded spectra of only moderate S IN, are rejec t ed. As 
in methods 1-3, damping constants from Barkl enTet al.l (j2000) 
were used when availabl e. Abundanc es were calculated with the 
MOOG 2002 package dSnedenll 19731 ). 

The effective temperature is determined by requiring excita- 
tion equilibrium for the set of Fe i lines with measured equiva- 
lent widths less than 130 mA. The microturbulence velocity (£ t ) 
is solved for in the standard way. As a first pass, models with sur- 
face gravity of log g = 4.5, and with solar metallicity are used. 
The [Fe/H] value is then determined from the set of Fe i lines 
with lower excitation potential greater than 4.0 eV and equiv- 
alent widths less than 130 mA using that choice for T e ff. The 
calculation is repeated with different metallicity models if the 
output of the first pass differs significantly from that of the input 
model atmosphere. The surface gravity then follows by demand- 
ing ionisation equilibrium between Fe i and Fe n. The rationale 
behind this method is further described in Cohen et al. (2009, in 
prepar ation), and additional details can be found in Cohe n~et"aT1 
d2008h . 

3.6. Comparisons of stellar parameters from methods 1-4 

Below we present the stellar parameters for the four stars and 
compare them to the ones found b y the different met hods and by 
iJohnson et al] (I2007L l2008h and lCohen etail d2008). All values 
are given in Table [3] 



4 The interpolator by T. Masseron, and the new MARCS models are 
available on the MARCS website http : / /marcs . astro .uu.se 
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Table 2. Bensby's linelist with measured equivalent widths and calculated elemental abundances for OGLE-2008-BLG-209S. Also 
included are the equivalent widths and elemental abundances measured by Bensby in the solar spectrum and in spectra for OGLE- 
2006-BLG-265S, MOA-2006-BLG-099S, and OGLE-2007-BLG-349S provided by Cohen and Johnson. For each line we give the 
\oggf value, lower excitation potential (x\), measured equivalent widths (Wa), derived absolute abundance (loge(X)), abundance 
normalised to the Sun (on a line-by-line basis) ([X/H]). The "flag" for the Sun indicates that a line has not been measured in 
the solar spectrum and that the abundance given (and used in the normalisation) for that line is the average abundance based on 
all the other lines, that could be measured in the solar spectrum, of the same species. The table is only available in the online 
version of the paper and in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (138.79.125.5) or via 
http : //cdsweb . u-strasbg . fr/ Abstract . html. 



| Sun | OGLE-2008-BLG-209 | OGLE-2006-BLG-265 | MOA-2006-BLG-099 | OGLE-2007-BLG-349 

Element Wavelength logy/ X \ | loge(X) a flag] W A loge(X) [X/H] \ W A loge(X) [X/H] ] W A loge(X) [X/H] \ W, t loge(X) [X/H] 



Al 1 


5557.063 


-2.21 


3.14 


10.4 


6.44 


17.1 


6.45 


0.01 


39.9 


7.04 


0.60 


29.4 


7.00 


0.56 


37.9 6.89 


0.45 


Al 1 


6696.023 


-1.62 


3.14 


44.9 


6.62 


47.9 


6.45 


-0.17 


87.8 


7.09 


0.47 


58.8 


6.84 


0.23 


89.5 7.03 


0.41 



Table 3. Stellar parameters as determined from the different methods. 



Object 


(V-I)o 


/o 


Vo 


My 


d 


Teg 


log g 


it 


Mass 


log 


e(Fe) 


[Fe/H] 


Age 






[mag] 


[mag] 


[mag] 


[mag] 


[kpc] 


[K] 


[cgs] 


[kins-'] 


[M Q ] 


Fei 


Fen 




[Gyr] 




OGLE-2008-BLG-209S 






18.35 


3.83 


8.0 


5243 


3.82 


1.01 


1.02 


7.24 


7.24 


-0.33 


8.5 


Method 1 




0.73 


17.05 


17.68 


3.16 


8.0 


5307 


3.69 


1.13 


1.31 


7.26 


7.12 


-0.31 


3.8 


Method 2 














5250 


3.60 


1.40 




7.04 


7.04 


-0.42 




Method 3 














5325 


3.60 


1.30 




7.16 


7.16 


-0.32 




Method 4 


OGLE-2006-BLG-265S 






19.10 


4.58 


8.0 


5486 


4.24 


1.17 


1.05 


8.02 


8.01 


0.44 


7.5 


Method 1 




0.68 


18.11 


18.69 


4.17 


8.0 


5526 


4.11 


1.25 


1.10 


8.02 


7.91 


0.44 


7.0 


Method 2 










4.30 




5650 


4.40 


1.20 




8.05 


8.07 


0.56 




Johnson et al. (2007) 


MOA-2006-BLG-099S 






19.35 


4.83 


8.0 


5741 


4.47 


0.84 


1.12 


7.96 


7.95 


0.39 


* 


Method 1 




0.74 


18.17 


18.81 


4.29 


8.0 


5852 


4.32 


1.03 


1.18 


7.99 


7.81 


0.42 


2.8 


Method 2 










4.50 




5800 


4.40 


1.50 








0.36 




Johnson et al. (20081 


OGLE-2007-BLG-349S 






19.35 


4.83 


8.0 


5229 


4.18 


0.78 


0.94 


7.96 


7.95 


0.41 


13.0 


Method 1 




0.78 


18.72 


19.40 


4.88 


8.0 


5210 


4.18 


0.76 


0.93 


7.96 


7.96 


0.40 


13.7 


Method 2 














5480 


4.50 


1.00 




8.02 




0.51 




Cohen et al. (20081 



Table 4. Elemental abundance ratios, [X/H], as determined from the different methods. 



Object Oi Nai Mgi Ah Sii Cai Til Tin Cri Cm Nil Fei Fen Zm Yn Ban 



OGLE-2008-BLG-209S 0.06 -0.18 0.01 -0.04 -0.13 -0.15 -0.14 0.03 -0.29 -0.27 -0.28 -0.33 -0.28 -0.31 -0.45 -0.30 Method 1 
-0.07 -0.12 0.09 0.00 -0.15 -0.08 -0.09 -0.05 -0.23 -0.37 -0.27 -0.31 -0.40 -0.37 -0.54 -0.32 Method 2 
-0.08 -0.28 -0.14 -0.09 -0.25 -0.30 -0.26 -0.04 -0.33 - -0.36 -0.42 -0.44 -0.36 -0.77 Method 3 

-0.32 Method 4 



OGLE-2006-BLG-265S 


0.28 
0.18 
-0.03 


0.58 
0.64 
0.62 


0.56 
0.63 
0.48 


0.52 
0.55 
0.81 


0.49 
0.47 
0.48 


0.37 
0.43 
0.33 


0.42 
0.45 
0.46 


0.34 
0.27 


0.45 
0.47 
0.64 


0.26 
0.18 


0.49 
0.48 
0.62 


0.44 
0.44 
0.55 


0.43 
0.33 
0.57 


0.32 
0.30 


0.52 
0.45 


0.29 Method 1 
0.25 Method 2 
0.53 Johnson et al. (2007) 


MOA-2006-BLG-099S 


0.37 


0.46 


0.45 


0.42 


0.37 


0.31 


0.36 


0.45 


0.25 


0.28 


0.39 


0.39 


0.41 


0.47 


0.42 


0.37 Method 1 




0.22 


0.56 


0.56 


0.48 


0.38 


0.42 


0.43 


0.36 


0.32 


0.16 


0.42 


0.42 


0.28 


0.44 


0.31 


0.36 Method 2 




0.20 


0.45 


0.52 


0.47 


0.42 


0.25 


0.31 


0.14 






0.38 


0.36 


0.18 


0.08 




-0.25 Johnson et al. (2008) 


OGLE-2007-BLG-349S 


0.34 
0.37 
0.14 


0.55 
0.54 
0.77 


0.43 
0.42 
0.59 


0.49 
0.48 
0.49 


0.47 
0.48 
0.47 


0.31 
0.30 
0.55 


0.43 
0.41 
0.53 


0.41 
0.42 
0.67 


0.40 
0.39 
0.53 


0.41 
0.42 


0.48 
0.48 
0.67 


0.41 
0.40 
0.51 


0.42 
0.44 
0.43 


0.49 
0.50 
0.42 


0.26 
0.27 


0.19 Method 1 
0.19 Method 2 
0.38 Cohen et al. (2008) 



OGLE-2008-BLG-209S: The derived stellar parameters, 
\ogg = 3.6 to 3.8, and T e $ = 5200 to 5300 K, are consistent 
with OGLE-2008-BLG-209S being a subgiant star which also 
can be seen in Fig. [8] Furthermore, we find that OGLE-2008- 
BLG-209S has a sub-solar metallicity of [Fe/H] * -0.33 
(ranging between [Fe/H] = -0.42 and [Fe/H] = -0.32 for the 
different methods). This relatively low metallicity is in stark 
contrast to the previous three Bulge dwarf stars that were all 
determined to have highly super-solar metallicities: two with 



[Fe/H] > 0.5 dJohnson et alj|2007t ICohen et al. 2008), and one 
with [Fe/H] > 0.3 dJohnson et al.fl2008l) . From methods 1 and 
2 we find an age for OGLE-2008-BLG-209S ranging between 
4-8 Gyr (see Fig.©. 

OGLE-2006-BLG-265S: Th e effective temperatu re from the 
different methods and from Joh nson et alj (120071) agree quite 
well and range between 5500 to 5650 K. Also the surface grav- 
ity is well constrained between \ogg = 4.1 to 4.4, values that 
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are typical for dwarf stars. Using methods 1 and 2 we find that 
OGLE-2006-BLG-265S is a very metal- rich star at [Fe/H] = 
0.44. It is however 0.12dex lower than in lJohnson et alj (2007) 
who finds [Fe/H] = 0.56. From methods 1 and 2 we find that 
OGLE-2006-BLG-265S has an age of approximately 7 Gyr (see 
Fig. H. From Fig. we see that OGLE-2006-BLG-265S is a 
dwarf star close to the main sequence turn-off. 



MOA-2006-BLG-099S: Methods 1 and 2 give r eff = 5741 to 
5852 K, log g = 4.32 to 4.47, and [Fe/H] = 0.39 to 0.42, which 
are va lues that are very similar to the results by iJohnson et al.1 
(120081) that found r eff = 5800 K and [Fe/H] = 0.36. While 
Johnson et al. (2008) used the microlensing technique to esti- 
mate My (and thus log g), they determined the effective temper- 
ature from profile fitting of the wings of the Ha and H/3 Balmer 
lines. Since MOA-2006-BLG-099S falls just outside the limits 
of the isochrones we can only give an upper limit to its age, 
which should be approximately 3 Gyr (see Fig. [8j. Given the 
stellar parameters and the position in the My-T e f{ diagram, it is 
clear that MOA-2006-BLG-099S is a main sequence dwarf star. 



OGLE-2007-BLG-349S: The stellar parameters based on meth- 
ods 1 and 2 are very similar: T e ff sa 5220 K and [F e/H] » 0.40. 
These results are lower than the ones found by ICohen et al.l 
d2008h . r eff by ~ 260 K and [Fe/Hl by ~ .1 dex. Given the un- 
certainties claimed by ICohen et al.l d2008l) of 100K in r e fj and 
~ 0.1 dex in [Fe/H], and the uncertainties for methods 1 and 2 as 
listed in Sect. [5] we see that the metallicities are compatible and 
also the temperatures, but just barely. We find that OGLE-2007- 
BLG-349S is a dwarf star close to the main sequence turn-off, 
and that it has a high age of approximately 13 Gyr (see Fig . [8} . 



In summary, the differences we find in the stellar parameters 
among the groups (and methods), with the possible exception 
of OGLE-2007-BLG-349S, are easily within the error bars (see 
Sect. |5]l. 



3.7. Inconsistencies 

Equivalent width measurements and log gf -values: For OGLE- 
2008-BLG-209S, method 3 gives a T e g similar to methods 1 and 
2 but a lower metallicity: 0.2 dex in absolute values and 0. 1 dex 
when normalised to the Sun. The decrease in the difference when 
normalising to the Sun can be explained as a result of the lower 
Solar abundances that method 3 produces compared to methods 
1 and 2 (see Sect.|4]and Fig.|6]l. The difference of 0.2 dex in the 
absolute Fe abundances can be explained by difference of 0.2 dex 
in log g and 0.4 km s" 1 in £. In lBensbv etall d2003ll2005h where 
we in detail investigated the effects of uncertainties in these pa- 
rameters we see that the above differences in £ t and log g would 
result in a changes of ~ 0.10 - 0.15 dex and ~ 0.03 dex, respec- 
tively, in the output Fe abundance. Therefore it is clear that it 
is mainly the microturbulence parameter that is to blame for the 
different absolute Fe abundances between methods 1 and 3. In 
Sect. |4] we also show that the output abundances from methods 1 
and 3 are very similar if the input equivalent widths are similar 
and the model stellar parameters are the same. Hence, there are 
no problems with the model stellar atmospheres (MARCS 1997 
versus MARCS 2008), nor the different abundance programs 
(EQWIDTH versus TurboSpectrum), that are used by methods 
1 and 3 in the abundance analysis. 



Method 4 gives a higher r eff for OGLE-2008-BLG-209S 
than methods 1, 2, and 3, and a similar metallicity. As shown 
in Sect. 13. ll there are significant differences between the Fe i and 
Fe ii equivalent widths as measured by Bensby and Cohen. The 
Bensby measurements are 4 % and 12 % larger for Fe i and Fe n 
lines, respectively, and can possibly affect the accuracy of the 
ionisation balance that are used in finding the stellar parameters. 
To investigate if this could explain the differences in the stellar 
parameters between methods 1 and 4 we divide the Fe i and Fe n 
equivalent widths of Bensby by 1.04 and 1.12, respectively, and 
redo the determination of the stellar parameters. The effect is a 
slightly lower effective temperature, 5214K, from method 1 and 
does not resolve the problem with the different effective temper- 
atures between methods 1 and 4. However, they only differ by 
80 K, which is well within the uncertainties. 



In Sect. |3.1| we also found differences in the log gf- values for 
the Fei and Fen lines between the Bensby and Cohen linelists. 
Average differences are -0.03 + 0.07 dex for F ei lines and 
-0.05 ±0.09 dex for Fe n lines, with the values from lCohen et al.l 
(2008) being the larger ones. In methods 1 and 2 the abundances 
from Fe i and Fe n lines have been used in their absolute form, 
i.e. they have not been normalised to the Sun, when determining 
the stellar parameters. Therefore, if the different methods/studies 
have different log ^/-values for the Fe i and Fe n lines it may 
have an impact on the stellar parameters. Applying similar cor- 
rections to the Fe i and Fe n log ^/-values of the Bensby linelist, 
as was done above to the Bensby equivalent widths, and redoing 
the analysis have essentially no effect on the derived parameters. 

Furthermore, the recent study by Melendez & B arbuvl(l2009l) 
cl aims that there might be pr oblems with the Fe n log g /-values 
oflRaassen & Uvlingsl {T998) that were adopted by Bensbv ^t al.l 
(12001" and that we use i n methods 1 and 2. According to 
iMelendez & Barbuvl d2009t) these log ^/-values are both inaccu- 
rate and imprecise. We therefore check our stellar parameters for 
OGLE-2008-BLG-209S us ing the revised Fe n log gf- values by 
IMelendez & Barbuvl (120091) and we find slightly different values 
for the stellar parameters: = 5 197 K, log g = 3.88, and £ t = 
0.93. The age of the star also becomes higher, 12.5 + 6 Gyr com- 
pared to 8.5+6 Gyr. Apart from the age difference, these changes 
are really marginal. Therefore , and a s about half the loggf- 
values by Me lendez & Barbuvl (|2009) are based on an inverse 
solar analysis (which most lik ely is the cause for the cla imed 
abundance spread decrease by Melendez & Barbuy 2009 0, w e 
will keep the theoretical values by iRaassen & Uvlings (1998). 
We will then have a single source fo r our log g f- values, and wil l 
also be independent on the methods M elendez & Barbuvl {2009) 
use to measure equivalent widths for these lines, and to their 
choice of model stellar atmospheres used in the inverse solar 
analysis. 



Colou rs and temperatures: From the calibration by lAlonso et al.l 
d 19961) we check what temperature we should expect given the 
inferred colour of OGLE-2008-BLG-209S. The Alonso cali- 
bration is however in the Johnson system only, so we con- 
vert the Johnson-C ousin co l our of OGLE-2008-BLG-209S us- 
ing the relation by iFernid d 1983b giving (V - T)j = 0.924. 
Applying Eq. (6) in lAlonso et al.1 d 19961) we then get a tem- 
perature of 5564 K. And what colour would be predicted for 
this star at T=5250? Using the same equations (in reverse) we 
get a Johnson-Cousin colour (V — 7)q = 0.83. Double-checking 
with the colour- r e ff calibration of Ramirez & Melendezl ([2005) 
we find that this star should have an intrinsic Johnson-Cousin 
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Fig. 6. Differences in measured equivalent widths and elemen- 
tal abundances for the Sun for the Bensby (B) and Johnson (J) 
analyses of the Ganymede solar spectrum, a) shows the On av- 
erage the Bensby equivalent widths are 0.4 ±10% larger than 
Johnson's (see plot a), and the corresponding Bensby abun- 
dances are on average 0.077 ± 0.095 dex lower than Johnson's 
(see plot b). Plot c) then shows the differences when the Bensby 
abundances are computed using exactly the same stellar param- 
eters as used by Johnson (see discussion in text). The differences 
then decrease to less than 0.01 dex. Differences and standard 
deviations (in parentheses) are indicated in the figures and are 
based on 128 spectral lines that Bensby and Johnson have in 
common. 

(V - /)o ~ 0.82 if its temperature is about 5250. Hence it seems 
that all is fine and the temperatures we derive should be good. 

However, an intrinsic colour (V-/)o ~ 0.82 means that there 
is an inconsistency in method 2, because T e ff is adopted in dis- 
agreement of this colour (0.73 compared to 0.82). This could 
indicate that the reddening perhaps is wrong. 

But, is this discrepancy possible? In principle yes. We do 
have an example of a source (OGLE-2008-BLG-513S) that has 
0.3 mag less extinction than the clump. So 0.1 mag is not un- 
heard of. Nevertheless, this discrepancy is still bigger than what 
has been seen for the other three dwarf stars. Anyway, given the 
estimated uncertainties in r e jf (see Sect. [5} it seems that an effec- 
tive temperature of 5243 K is not unreasonable for OGLE-2008- 
BLG-209S. 

4. Solar analysis 

For methods 1 and 2 a solar analysis was performed on a 
spectrum of Jupiter's moon Ganymede that was obtained in 
March 2007 with the MIKE spectrograph using the same in- 
strument settings as for the observations of OGLE-2008-BLG- 
209S. Analysing the solar spectrum in a similar way (as we know 
the absolute magnitude of the Sun we do not require ionisation 
equilibrium) as for methods 1 and 2, we derive T e s = 5790 K, 
loggr = 4.45, £ = 1.04 kms" 1 , and log e(Fe) = 7.56 for the 
Sun. The final elemental abundances based on methods 1 and 2, 
will be normalised to those of the Sun from this analysis. The 



normalisation will be done on a line-by-line basis, and then av- 
eraged for each element, making the results strictly differential 
to the Sun. This way of normalising the abundances neutralises 
uncertainties and errors in the log gf- values. It should be noted 
that when a line in the Solar spectrum could not be measured, 
or, in the case of Fe lines, the line strength exceeded 90 mA, the 
average abundance from all the other lines of the same species 
were used for the solar abundance of that line. These cases are 
marked by "flag=l" in Table|2] 

In Table [5] we give the average solar abundances for each 
species as derived by Bensby from the MIKE Ganymede spec- 
trum. For comparison purposes we al s o give the standard Solar 
abund ances as given by Aspl und et al.l (120051) and lGrevesse et al] 
(120071) . Generally our Solar abundances agree within 0. 1 dex of 
the standard ones, with O and Ba being the exceptions. The dif- 
ferences, all though they are small, illustrates the importance 
of doing a strictly differential analysis toward the Sun, deriv- 
ing your own Solar abundances instead of adopting the tabulated 
ones. 

For method 3, a solar analysis was performed using the same 
spectrum of Ganymede as was used for methods 1 and 2. The 
line list and line parameters were the same, and a model stel- 
lar atmosphere for the Sun was interpolated in the same grid of 
model stellar atmospheres as was used for OGLE-2008-BLG- 
209S, using r eff = 5770, logg = 4.44, [Fe/H] = (log e(Fe) = 
7.45), and £ = 1.1 kms" 1 . 

For method 4, the same set of lines as used for OGLE- 
2008-BLG-209S we re analysed in the Solar spectrutr0 by 
IWallace et al.l (119981) . These were then use to determine [X/H] 
for each species for OGLE-2008-BLG-209S. 

As the same Ganymede spectrum was used by both Bensby 
and Johnson to determine the Solar elemental abundances we 
show in Figure U a comparison of measured equivalent widths 
and absolute elemental abundances for 128 spectral lines in com- 
mon. The equivalent width measurements are in good agreement 
and differs by only 0.5%, with the measurements by Bensby be- 
ing the slightly larger ones. Surprisingly, this good agreement 
in the equivalent widths is not reflected in the derived elemen- 
tal abundances for the Sun. On average the abundances from 
Bensby's analysis are 0.08 dex lower than the abundances from 
Johnson's analysis (see Fig. [(J)). This is puzzling, but it turns out 
that the solar models used by Bensby and Johnson differ. While 
Johnson used a MARCS model with log e(Fe) = 7.45, Bensby 
used a MARCS model with log e(Fe) = 7.56, and of course the 
values for the microturbulence parameter also differed slightly. 
So, redoing the Bensby solar abundances, using exactly the same 
set of model stellar parameters as used by Johnson, the differ- 
ences in the solar abundances become less than 0.01 dex (see 
FigS). This also explains why the differences between methods 
1 and 3 in the Fe abundances are different depending on if com- 
parisons are made between absolute abundances or normalised 
abundances (see Table|3]l. As we are primarily interested in com- 
paring stars with each other it is [] -notation that is relevant here 
and it is encouraging to see how well the different methods do 
reproduce the results (once any differences in the methodology 
are fully accounted for). This is perhaps not surprising but worth 
noting and serves as a reminder that if one attempts to combine 
data from different studies then careful normalisation is a must. 

Equivalent widths and absolute elemental abundances (also 
for the Sun) and normalised abundances are given for individual 
lines are given in Table [2] Average elemental abundance ratios 
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Table 5. Our solar abundances based on observations of 
Jupiter's moon Ganymede. Last c olumn give the standa rd Solar 
photospheric abundances given bv lAsplund et al.l d2005l) . 



Ion 


log £(X) 


N 


A2005 


Fe i 


7.57 ± 0.08 


148 


7.45 


Fe ii 


7.52 ± 0.13 


17 


7.45 


Oi 


8.86 ±0.02 


3 


8.66 


Nai 


6.29 ± 0.08 


4 


6.17 


Mgi 


7.62 ± 0.08 


5 


7.53 


All 


6.48 ± 0.07 


7 


6.37 


Si i 


7.60 ± 0.07 


27 


7.51 


Cai 


6.32 ±0.11 


19 


6.31 


Til 


4.93 ±0.11 


22 


4.90 


Tin 


4.90 ±0.17 


14 


4.90 


Cri 


5.55 ± 0.06 


6 


5.64 


Cm 


5.69 ± 0.04 


5 


5.64 


Nil 


6.22 ± 0.07 


42 


6.23 


Zni 


4.60 ± 0.08 


3 


4.60 


Yn 


2.12 ± 0.1 1 


4 


2.21 


Ban 


2.42 ± 0.07 


4 


2.17 



are given in Table [4] Absolute and normalised abundances from 
the Fe i and Fe n lines are also given in Table [3] 
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5. Errors 

5.1. Random errors in stellar parameters 

We are performing standard LTE abundance analyses and the 
associated errors and uncertainties for the stellar parameters 
should be of standard natur e. Methods 1 and 2 are similar to 
the analysis carried out in iBensbv et al.l (120031 120051) where 
errors of 80 K in T e g and 0. 1 dex in log g are quoted. The 
line-to-line scatter in those studies for [Fe/H] are of the order 
0.08 dex. For OGLE-2008-BLG-209S the standard deviation of 
the mean abundance in the Fe i abundances are about double that, 
<x[Fe/H] = 0.14 (see Table [6]), which most likely is an effect of 
th e lower quality of the spectr um as compared to the spectra used 
in IBensbv et aU d2003Ll2005l) that had S /N > 200. Therefore we 
estimate that the uncertainties in T e ff and log g should be higher 
for OGLE-2008-BLG-209S. A conservative estimate is 0.2 dex 
in log g and 200 K in T e g. The same numbers should hold for 
OGLE-2006-BLG-265S and MOA-2006-BLG-099S that have 
errors in the Fei abundances comparable to OGLE-2008-BLG- 
209S (see Tabled). The spectrum of OGLE-2007-BLG-349S is 
however of significantly higher quality than the others, which 
also is reflected in the lower Fe i standard deviation around the 
mean abundance, 0.10 dex (see Table [6j. Hence, for OGLE- 
2007-BLG-349S we estimate that the errors in T e ff, log g are of 
the order 150 K and 0.15 dex, respectively. The magnitudes of 
these e rrors a re in line with the value s given by iJohnson et al.l 
J2007L 120081) and ICohen et al.l (|2008) for OGLE-2006-BLG- 
265S, MOA-2006-BLG-099S, and OGLE-2007-BLG-349S. 

In the last column of Table|6]we list by how much the [X/Fe] 
ratios are affected by these uncertainties in the stellar param- 
eters. The values a re t aken as double the values as listed in 
IBensbv et al.l d2005) and lBensbv et al.ld2004 (for [O/Fe]) where 
we investigated this in detail (see discussion above). 

5.2. Elemental abundances and the robustness of 
abundance ratios 

We find that the [X/Fe] ratios generally agree well between 
the different methods/studies, and it is essentially only oxygen, 



Fig. 7. Comparisons of derived [X/Fe] ratio s for method 1 (open 
circle s), m ethod 2 (filled circle s), and from l Johnson et all (2007. 
2008) and lCohen etail d2008l) (stars). F or OGLE-2008-BLG- 
209S the stars show the values from method 3. For oxygen two 
values are given side by side: on the left hand side the LTE value 
as it comes of the analysis, and on the right hand side the NLTE 
value corrected by t he empirical NLTE correction formula from 
IBensbv et al.ld2004l) . 



zinc, yttrium, and barium that show any real discrepancies for 
one or two sta rs. The big d i fferenc e in [O/Fe] that can be seen 
compared to IJohnson et aU d2007l) for OGLE-2006-BLG-265S 
is clearly an effect of the almost 200 K difference in effective 
temperature. 

The Ba abundances from IJohnson et al.l (|2007, 2008]) are 
based on line synthesis, and hence the equivalent widths can not 
be directly com pared. However, exc ept for MOA-2006-BLG- 
099S, where the Johnso net al.l d2008l) [Ba/Fe] ratio is extremely 
low, the agreement is good. The low value of MOA-2006-BLG- 
099S is also not recognised in OGLE -2006-BLG-265S, a s imi- 
larly metal-rich star. As mentioned in IJohnson et alj d2008l) the 
reason for the low [Ba/Fe] value of MOA-2006-BLG-099S was 
the high value they measured in the Sun. The good agreement 
for the other stars between Ba as derived by line synthesis and 
through equiva lent width measurem ent analysis illustrates the 
point raised by iBensbv et al.l d2005l) that because the dominant 
isotopes for Ba are the even isotopes, which do not have hy- 
perfine structure, the inclusion of HFS does not substantially 
change the abundances. Furthermore, the line synthesis of the 
Ban line at 6141 A done by Johnson in method 3 for OGLE- 
2008-BLG-209S included an Fei line that blends with the fea- 
ture. However, this line contributes very little to the line pro- 
file. A synthesis of the line showed that the equivalent width in- 
creased from 127.9 nm to 132.7 nm if the Fe 1 line was included. 
Therefore, using the equivalent width and assuming that it is all 
Ba leads to a very small change in the Ba abundance. 
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6. Bulge membership for OGLE-2008-BLG-209S 

Bulge microlensing sources tend to be slightly biased to be more 
distant th an the mean d i stance to the Bulge ("the Bulge clump"), 
see e.g. iKane & Sah u (2006). The lens must obviously be in 
front of the source, so if the lens is in the Bulge (which most 
are), the source must be drawn from the Bulge stars behind it, so 
on average more distant than the Bulge centre. In fact, the bias 
is stronger than that because the probability of lensing goes ba- 
sically as y/(d source - duns) for lenses and source in the Bulge. 
For disc lenses (the minority) there is no significant bias. This 
bias would tend to make the source more luminous than the 
microlens-model based estimate based on assuming it is at the 
same distance as the red clump. This goes in the opposite direc- 
tion of the difference between the spectroscopic and microlens 
determinations of My . 

Based on our spectroscopic parameters from method 3, 
r e ff = 5250 and \ogg = 3.6, we use the Y 2 -isochrones at 
Z = 0.007 to estimate an Mi of 2.5. Combined with the Iq 
1 7.05 microlensing estimate, we derive a distance to OGLE- 
2008-BLG-209S of 8kpc. Now, if we take My = 3.83 from 
Table[3] (method 1) and V = 17.68 (microlensing estimate from 
method 2), the distance comes out to 6kpc. The difference is 
probably due to the higher gravity (i.e., fainter absolute mag- 
nitude). Since the errors on the gravity are not insignificant (be- 
tween the small number of Fe n lines and the correlation between 
temperature and gravity, not to mention systematic uncertainties 
between isochrone and spectroscopic determinations), there are 
large errors attached to both these distances. 

However, the high heliocentric radial velocity of OGLE- 
2008-BLG-209S (see Sect. |2T4j» demonstrates that it is likely to 
be in the Bulge, and the relatively high radial velocities of the 
other three dwarf stars also strongly suggest Bulge membership. 
This high velocity dispersion of our four stars is also in agree- 
ment with what is seen in large samples of gia nt stars in the 
Bulge (e.g. lSadler et alJll996tlZoccali et al.ll2008l) . 

In conclusion, based on the radial velocity, probability of 
microlensing, and small distance from the Galactic plane, we 
expect this star to be in the Bulge. This is consistent with the 
range of distances derived from spectroscopic parameters and 
isochrones. 

7. Stellar age and stellar mass 

When the final values of \ogg, T e e, and £ t have been found 
we use the fundamental astrophysical relation g cc M. ■ T S /L 
to estimate the luminosity (and hence the apparent magni- 
tude), assuming a distance of 8kpc to the star. For methods 
1 and 2 stellar ages were then estimated using the Yonsei- 
Yale i sochrones (lYi et alj2001l:lKim et al.l2002tlDemarque et al.l 
120041) . Individual sets of isochrones with appropriate metallic- 
ities and a-enhancements (in the case of OGLE-2008-BLG- 
209S) were calculated for the different stars. Ages were then 
read off from the best fitting isochrone in the log r e ff-My plane 
(see Fig. |8). Upper and lower limits to the ages were estimated 
from the error bars based on the uncertainties in T e ff and My. 
For r eff we adopt error bars of 200 K for OGLE-2008-BLG- 
209S, OGLE-2006-BLG-265S, andMOA-2006-BLG-099S, and 
150 K for OGLE-2007-BLG-349S (see Sect. 0. The errors in 
My are difficult to estimate. Given the range of values by the 
different methods in Table [3] we adopt a value of 0.7 mag for the 
error in My for all four stars. This value is consistent with an 
uncertainty of 2kpc in the adopted distance (see Sect. |6]). These 
error bars are shown in the isochrone plots (see Fig. |8). 



Table 6. Final elemental abundance ratios, [X/Fe], based on 
method 1 . For each abundance ratio we also give the standard de- 
viation of the mean abundance (added in quadrature for the two 
elements used to construct the abundance ratio), and the number 
of lines used. These are the values used in subsequent abundance 
plots. The last column gives an estimate how the abundance ra- 
tios are affected by random errors in the stellar parameters. The 
adopted uncertainti es for the Bulge dw arf s tars are two times th e 
values as given in lBensbv et~all(l2005l) and lBensbv et all (12004). 
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[Y/Fe] 
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0.03 


-0.15 


-0.02 


-0.22 


0.12 




0.16 


0.17 


0.17 


0.17 






4 


3 


4 


4 





Also, stellar masses (AT) w ere determined from the evolu- 
tionary tracks bv lYi et aTJ d2003l) and are given in Table [3] 

8. Stellar populations 

Even if there are some discrepancies in the stellar parameters 
from the different methods that we have presented, they do pro- 
duce results that are very similar in terms of abundance ratios. 
This is an encouraging result as it is the abundance ratios that are 
the main results used in discussions of the enrichment of heavy 
elements in a stellar population. 

However, in order to put the new results for these four Bulge 
stars into context we will adopt the results from method 1. 
Method 1 is the method that most closely resembles the method 
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Fig. 8. The Yonsei-Yale (Y 2 ) isochrones (lYi et al.ll200"Tl iKim et all 120021; iDemarque et al.l 12004b that were used to estimate the 
stellar ages. Isochrones are plotted in 1 Gyr steps, with every 5 Gyr isochrone in dashed lines. For each target we used isochrones 
metallicities according to Table [3] and for OGLE-2008-BLG-209S also with an appropriate the a-enhancement. Filled circles 
represent the stellar parameters derived from method 1 and open circles the stellar parameters from method 2. Error bars indicate 
uncertainties of 0.7 mag in My and 200 K in T e ff. 



used in lBensbv etaT] d2003l [2005b and Bensby et al. (2009, in 
prep), and will enable a differential comparison of the Bulge 
stars with a large sample of dwarf and subgiant stars in the 
nearby thin and thick discs. This choice is for consistency, and 
importantly, it does not make any judgement regarding which set 
of abundances is the most accurate. 

8.1. Elemental abundance trends in the Bulge 

The Milky Way contains four major stellar populations; the halo, 
the thin and thick discs, and the bulge. 

In Fig. [9] we show the thin an d thick disc abundance trends 
from dBensbv et al.l I2003L I2005I) . The stars of the thick disc 
have higher a-to-iron abundance ratios than the stars of the 
thin disc. Specifically, the stars of the thick disc show a con- 
stant high [a/Fe] ratio for low metallicities and all the way up 
to [Fe/H] =s -0.4 after which a down-turn in [a/Fe] occurs, 
steadily declining towards solar values at solar metallicity. The 
stars of the thin disc, on the other hand, show a shallow decline 
in [a/Fe] at lower metallicities that levels out at solar values at 
solar metallicities, where they merge with the thick disc trends. 

In Fig.[9]we also show the results for the four dwarf/subgiant 
stars in the Bulge overplotted on the th i n and thick disc abun- 
dance trends from iBensbv et all d2003l I2005I) . Generally, the 
abundance ratios for OGLE-2008-BLG-209S are similar to what 
is seen for the thick disc at the same [Fe/H]. Interestingly, the 
metallicity of OGLE-2008-BLG-209S (-0.33 dex) is where the 
separation between the thin and thick discs is the greatest. If 
anything, this demonstrates that the observed [a/Fe] ratios in 
OGLE-2008-BLG-209S make it unlikely for it to be a thin disc 
star. 

The three other Bulge dwarf stars, at super-solar metallici- 
ties, show abundance ratios for the a-elements in close agree- 
ment with the metal-rich thin disc, i.e., close to the solar values 
([a/Fe] * 0). Under the assumption that all four stars are all 
genuine members of the Bulge, this means that the abundance 
trends of the Bulge, somewhere in between these points, should 
show a significant decrease in the a-to-iron ratios, signalling the 
onset of chemical enrichment from low-mass stars. 

Figure [10] shows the abundance trends for O, Mg, Ti, and Si 
in the Bulge as traced b y K and M giant stars (JRich & Origlia 
200j iRich etal] l2007t iLecureur et ail l2007t Fulbright et all 
2007; Melende z et alj2008l) . with the results for the microlensed 
dwarf and subgiant stars overplotted. The main chemical char- 
acteristics of the Bulge using giant stars can be summarised as 



follows: (1) The stellar formation history and chemical enrich- 
ment has been very fast. This is reflected by high [a/Fe] ratios 
even at solar metallicities; (2) The [O/Fe] ratio starts to decline 
at roughly [Fe/H] m -0.3, a signature of the onset of chemical 
enrichment by low-mass stars. (3) The [ O/Fe] -[Fe/H] trends i n 
the Bulge and the thick disc are similar (Me lendez et alj 2008). 
(4) Abundance trends based on other a-elements such as Mg, Si, 
and Ti are similar to the oxygen trend, with the exception that 
the a-enhancements for these elements, albeit slowly declining, 
tend to stay at somewhat higher levels with me tallicity as com- 
pared to oxygen (see, e.g. iFulbright et al1l2007l) . Also, there are 
significantly fewer giant stars from these studies where abun- 
dances of Mg, and especially Ti and Si have been determined 
as compared to oxygen. (5) Generally the abundance ratios from 
the giant stars and the dwarf/subgiant stars seem to agree. 



In order for a down-turn to occur in the first place in the 
[O/Fe] -[Fe/H] plot the oxygen enrichment (by massive stars) 
needs to decrease. The [O/Fe] ratio will then decrease, rela- 
tive to Fe production, with [Fe/H] as a result of the continued 
Fe enrichment from low-mass stars. The maximum possible de- 
crease of the [O/Fe] trend will occur if oxygen enrichment is 
shut off completely. The [O/Fe] ratio will then decrease by the 
same amount with which [Fe/H] increases. In the [0/Fe]-[Fe/H] 
plot this correspond to a line with a slope equal to - 1 . For sim- 
plicity, lets assume that this is the case. In Fig. [10] we show 
this line, and it is clear that, in order for the [O/Fe] ratio path 
to be able to decline to the observed [O/Fe] ratio in the three 
metal-rich Bulge dwarf stars, the oxygen production needs to be 
abruptly shut off at a metallicity not higher than [Fe/H] m —0.1. 
It is of course unlikely that the oxygen enrichment should turn 
off "over-night", instead it should be a gradual decrease, i.e. the 
slope of the line should be less than what we have indicated in 
Fig. [10] If so, it seems likely that the turn-over in the [O/Fe]- 
[Fe/H] trend as traced by dwarf and subgiant stars happens at 
approximately the same metallicity as for the g iant stars, i.e. at 
[Fe/H] * -0.3 to -0.4 dMelendez et al J 120081 This also hap- 
pens to be wh ere we see the "knee" in t he thi ck disc [a/Fe]- 
[Fe/H] trends dBensbv etalj|2003l 12004 120051) . Similar lines 
have been drawn in the Mg, Si, and Ti plots in Fig. [10] indicating 
that the enrichment from low-mass stars started at [Fe/H] w 
(under the assumption that the enrichment of these elements 
from massive stars turned off abruptly). 
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Fig. 9. Abundance trends , [X/Fe] versus [Fe/H], for the thin and thick discs (taken from lBensbv et al . 2003, 2005). Thin and thick 
disc stars are marked by open and filled circles, respectively. OGLE-2008-BLG-209S is marked by a red circle, OGLE-2006-BLG- 
265S by a yellow circle, MOA-2006-BLG-099S by a blue circle, and OGLE-2007-BLG-349S by a green circle. The error bars have 
two cross bars. The inner one represent the standard deviation of the mean abundance, and the outer cross bar represent the total 
error, i.e. standard deviation of the mean abundance and the uncertainty due to random errors in the stellar parameters, as given in 
Table |6j added in quadrature. 
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Fig. 10. Elemental abundance trends in the Bulge based on giant stars (marked by crosses ) (O, Mg, Si, and Ti fr om lRich & O riglia 
2005i lRich et al ll2007l:lFulbright et al|l2007t O and Mg from lLecureur et al.ll2007t O from lMelendez et alj|2008l) . Nearby thick disc 
giant stars from Melend ez et al.l d2008l) are marked by filled circles. The four stars in this study are marked as in earlier figures. The 
diagonal line marks the maximum slope that the [X/Fe] ratio can decrease in order to reach the levels of the three metal-rich dwarf 
stars (the actual slope is however likely to be shallower, see discussion in Sect. 18.11 ) 



8.2. A comparison of Bulge and inner disc dwarf stars 



There is currently no in situ stellar sample that traces the in- 
ner disc regions of the Galaxy using detailed elemental abun- 
dances. A possible way to form a compatible sample would be 
to sample high-velocity stars with orbits that make them likely 
to have come from the inner parts of the Galactic disc. For 
instance, Pompeia et al. ( 2002) studied a sample of, so called, 
nearby "bulge-like" dwarf stars. These were defined as stars with 
highly eccentric orbits (<? > 0.25), that do not reach more than 
1 kpc from the Galactic plane (Z max < 1 kpc), and that have a 
maximum peri-galactic distance {R p ) of 3 to 4 kpc. In Bensby et 
al. (2009, in prep) we have ~ 700 stars and 82 of tho se stars ful- 
fil the "bulge-like" criteria by Pompei a"et al.l d2002l) . Figure [TT] 
shows the [Ti/Fe]-[Fe/H] trends for these stars together with the 
four microlensed Bulge dwarf/subgiant stars. The observed trend 
for the "bulge-like" disc sample is very well defined with very 
little scatter. At sub-solar [Fe/H] there are no signs of [Ti/Fe] 
ratios seen in the local thin disc. This might not come as an sur- 
prise as we are picking stars on highly eccentric orbits, and if 
we were to trace the thin disc trend closer to the Galactic cen- 
tre, those stars would likely be on on almost circular orbits and 
hence never cross the Solar orbit. It could also be that these disc 
stars on highly eccentric orbits belong to an inner disc popula- 



tion, distinct from the local thin and thick discs. This will be 
investigated further in Bensby et al. (2009, in prep.). 

Plotting the microlensed Bulge dwarf and subgiant stars on 
top of the "bulge-lik e" disc sample (using the definition from 
iPompeia et afllfool in Fig. E]we see that OGLE-2008-BLG- 
209S nicely fits into the trend and that the other three Bulge 
dwarf stars form a metal-rich extension of the trend. 



8.3. The metallicity distribution of the Bulge 

In Fig. [T2]we plot the metallic ity distribution of th e 521 giant 
stars in the Bulge observed by IZoccali et alj (12008) in Baade's 
window (l,b) = (1.14,-4.18), and two other Bulge fields at 
(l,b) = (0.21,-6.02) and (l,b) = (0.0,-12. 0). Also shown 
is the MDF of the 24 giant stars analysed by Fulbright et al. 
d2007l) . on which we have overplotted the metallicities of the 
four Bulg e stars analy s ed in t his work. Note that the distribution 
of the Fulbrig ht et alj d2007) is not representative of the MDF 
in the Bulge. Their sample is designed to probe the full range 
of [Fe/H], not to be a representative sample. As can be seen 
the Bulge MDF based on the giant stars range from very low 
metallicities, [Fe/H] w -1.3, to super-solar values comparable 
to the three previously published dwarf stars at [Fe/H] « 0.4- 
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Fig. 11. 82 nearby dwarf stars, selected from Bensby et 
al. (2009, in prep.), that have kinematic properties making them 
likely to come from the in ner Galactic disc/Bulge . The kinematic 
criteria are the same as in iPompeia et al.l d2002), i.e., e > 0.25, 
Zmax < 1.0 kpc, and R min < 4kpc. OGLE-2008-BLG-209S 
(red circle), OGLE-2006-BLG-265S (yellow circle), and MOA- 
2006-BLG-099S (blue circle) have also been included. 



0.5 dJohnson et alJl2007l 12001: ICohen et al.ll2008h . The average 
metallicity is [Fe/H] * -0.2 (IZoccali et alj|2008ll 

An issues that has been discussed against results based on 
dwarf stars that have been microlensed is that the microlensing 
is affecting the spectrum of the source star and therefore we are 
not getting the right answer. This is why these stars (for some 
reason) are all very metal-rich ( IZoccali et al.l 12.008b . OGLE- 
2008-BLG-209S, the most recently analysed microlensed star, 
although being a subgiant, has a metallicity of [Fe/H] = -0.33, 
and proves that it is indeed possible to find metal-poor mi- 
crolensed stars. Also, since OGLE-2008-BLG-209S is approx- 
imately at the distance of the stars in Baade's Window (because 
it's very close spatially to BW, see Fig. [U), its metallicity is very 
compatible with the MDF of giant stars if they are spatially co- 
incident (whic h is the case of OG LE-2008-BLG-209S and the 
giant stars by IZoccali et alj |2008). So this strengthens the idea 
that the analysis of microlensed dwarf and subgiant stars needs 
to be taken seriously. 

8.4. Formation and chemical history of the Bulge 

There are essentially three main formation scenarios for how 
galactic bulges may form. 

Firstly, we have the classical scenario where the bulge form 
in the final phases of the m onolithic collapse of the protogalactic 
cloud (Eggenet al. 1962). First the halo forms and as the col- 
lapse progresses also the bulge. The similarity between the dis- 
tributions of angular mom entum in the Bulge and halo suggest 
that this could be the case (TWyse & Gilmordl 19921) . Secondly, it 
is possible that bulges form from mergers in a ACDM Universe. 

The outcome from these two scenarios is that the abun- 
dance trends in the Bulge could show many types of signatures. 
Depending on the star formation history, it is possible to get high 
a-to-iron ratios extending to super-solar metallicities. However, 
as both of these scenarios happen in the very early phases of the 
formation history of the Galaxy, they would result in a Bulge 
that contain very old stars, as old as those found in the stellar 
halo. Even if not exclus ively so, old stars are generally what is 
found in the Bulge (e.g. Feltzing & Gilmorel2000tlOrtolani et al.l 
[19951) . However, in the collapse model, the average metallic- 
ity should increase as the collapse progresses inwards, resulting 
in metallicity gradients. Tentative evidenc e for a small metal - 
licity gradient in the Bulge was found bv IZoccali et ail f2008). 
see Fig. Q~2] For the hierarchical merger scenario in a ACDM 
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Fig. 12. Metallicity distribution of the Bulge. Top panel 
shows metalli c ity h istogram of the 24 giant stars from 
Fulbright et al.l (120071) . on which the metallicities of the four 
Bulge dwarf/subgiant stars analysed in this work, marked as in 
previous plots (arbitrarily shifted vertically). The three bottom 
panels sho w the MDF of the g iant stars in three different Bulge 
fields from lZoccali et all d2008l). The (I. b) = (1.14, -4.18) field 
has 244 stars, the (l,b) = (0.21,-6.02) field 213 stars, and the 
(l,b) = (0.0, -12.0) field 104 stars. 



Universe, where the in-falling "pieces" will be randomly dis- 
tributed within the Bulge, no metallicity gradients are expected. 

The third scenario is that bulges form through secular evolu- 
tion, in which a bulge-like structure is built up through the slow 
rearrangement of energy and mass as a result of interactions be- 
tween for instance a galactic bar and a galactic disc. If secular 
evolution dominates then it is quite likely that abundance trends 
in the bulge would mimic those of the central galactic disc as 
that is where the material and energy would come from to form 
the bulge. For the Milky Way bulge, the secular scenario has 
until recently been seen as quite unlikely because of the differ- 
ing abundance patte rns in the Bulge to those in the disc (e.g . 
iFulbright et~aT1l2007l) . However, recently iMelendez et all ((2008) 
carried out a truly differential study of Bulge and thin and thick 
disc giant stars and found that they are indeed very similar; the 
Bulge almost perfectly mimics the oxygen abundance trends of 
the nearby thick disc. Our results do not contradict this similar- 
ity. 

As also pointed out by Melen dez et al.l d2008l) . the simi- 
larity of the Bulge and the nearby thick disc oxygen trends, 
even though these two population clearly are spatially sepa- 
rated within the Galaxy, indicates that they have had very simi- 
lar, but not necessarily the same, chemical enrichment histories. 
Observations of dist ant galaxies suggest that the secular scenario 
is very common dGenzel et al.l 120081) . and with the new abun- 
dance results we find that it is also a plausible scenario for the 
Milky Way bulge. 
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9. Summary 

We have performed a detailed elemental abundance analy- 
sis of OGLE-2008-BLG-209S, a microlensed subgiant star in 
the Bulge. In contrast to the three previously published mi- 
crolensed Bulge dwarf stars, that all turned out to be extremely 
metal-rich, having metallic i ties greater t han [ Fe/H] > 0.35 
dJohnson et alJ 120071 120081: ICohen et alJ |2008), OGLE-2008- 
BLG-209S has a more modest sub-solar metallicity of [Fe/H] m 
-0.33. Interestingly, the metallicity of OGLE-2008-BLG-209S 
is where the separation in the [a/Fe] ratio between the thin and 
thick disc stars is the greatest. Furthermore, the abundance pat- 
tern of OGLE-2008-BLG-209S is similar to what is found in 
nearby thick disc dwarf stars and what is found in giant stars in 
the Bulge, i.e. enhanced a-to-iron ratios (e.g., O, Mg, Si, Ca, Ti). 

We have also carried out a re-analysis of the three mi- 
crolen s ed dw arf s tars, previou s ly pub lished by iJohnson et alJ 
d2007l 120081) and ICohen et alJ ((2008). This was done in or- 
der to get a homogeneous sample of Bulge dwarf and subgiant 
stars analysed with the same methods as for t he large sam- 
ple of nearby thin and thick disc dwarf stars by iBensbv et alJ 
(I2003L 120051) . With some exceptions we see good agreement 
in abundances between our results for MOA-2006-BLG-099S, 
OGLE-2006-BLG-265S, OGLE-2007-BLG-349 S and those in 
IJohnson et ail d2007l 120081) : ICohen etalJ d2008). For instance, 
the metallicity of OGLE-2006-BL G-265S is lower in ou r new 
analysis ([Fe/H] = +0.44) than what lJohnson et al.l d2007l) found 
([Fe/H] = +0.56). 

The three stars for which we could estimate ages are (us- 
ing method 1) 8.5 Gyr, 7.5 Gyr, and 13Gyr old, respectively. 
These ages are higher than what we see for thin disc stars at 
these metallicities. However, the stars of the thick disc are found 
to be older than the stars of the thin disc. This could suggest 
that the metal-rich stars of the Bulge have ages comparable to 
what we see in the metal-rich thick disc (at solar metallicities, 
IBensbv et al.ll2007b . It should be cautioned that individual ages 
based on isochrones could be subject to substantial systematic 
and random uncertainties. Therefore, ages for many more bulge 
dwarf and subgiant stars must be determined before firm conclu- 
sions can be drawn. 

We have also considered several methods in analysing the 
stars and find that elemental abundance ratios based on the meth- 
ods are in good agreement, showing small or negligible differ- 
ences. We would like to stress that the differences between the 
methods we used to get abundances from Bulge dwarf stars have 
been quantified, and that our investigation shows that one always 
should be very careful in combining information from differnet 
studies unless all "facts" are clearly explained (such as the solar 
abundances used for normalisation, etc.). 

The above results are based on less than a handful of 
dwarf/subgiant stars. It is very important to further confirm these 
results with more observations. However, high-magnification 
microlensing events, that currently seem to be the only feasible 
way of getting good quality high-resolution data of Bulge dwarf 
stars, are very rare. Out of the yearly thousands of detected mi- 
crolensing events from the OGLE monitoring of the Bulge, only 
a dozen are dwarf events. As it is impossible to know beforehand 
when these will occur, it is necessary to have a flexible observing 
schedule to be able to catch them. As such, the Rapid Response 
Mode with the UVES spectrograph on the Very Large Telescope 
should be optimal. 
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